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ABSTRACT 

Context. LS 5039/RX J1826. 2-1450 is one of the few High Mass X-ray binary systems from which radio and high energy TeV 
emission has been observed. Moreover, variability of the TeV emission with orbital period was detected. 

Aims. We investigate the hard X-ray (25 - 200 keV) spectral and timing properties of the source with the monitoring IBIS/ISGRI 
instrument on-board the INTEGRAL satellite. 

Methods. We present the analysis of INTEGRAL observations for a total of about 3 Msec exposure time, including both public data 
and data from the Key Programme. We search for flux and spectral variability related to the orbital phase. 

Results. The source is observed to emit from 25 up to 200 keV and the emission is concentrated around inferior conjunction. Orbital 
variability in the hard X-ray band is detected and established to be in phase with the orbitally modulated TeV emission observed with 
H.E.S.S. For this energy range we determine an average flux for the inferior conjunction phase interval of (3.54 ± 2.30) x 10"" erg 
cm"^ s"' , and a flux upper limit for the superior conjunction phase interval of 1.45 x 10"' ' erg cm"^ s"' (90% conf. level respectively). 
The spectrum for the inferior conjunction phase interval follows a power law with an index T = 2.0^^2 (90% conf. level). 
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1. Introduction 

; LS 5039/RX J1826.2-1450 is a High Mass X-ray Binary 
(HMXRB) system, which was discovered in the ROSAT all- 
sky survey_and identifi ed with optical follow-up observations 
, by iMotch et al.l ( Il997b . The compact companion is estimated 
■ to have a mass of M ~ 3.7 Mq, orbiting a bright, V = 11.2, 
06.5V((f)) star with an orbital period Poi-b ~ 3.9060 + 0.00017d, 
an eccentrici ty e - . 35 + 0.04 and a distance of c/ = 2.5 + 
0.1 kpc ( Casares et al.L l2005i) . Based on the observation of per- 
sistent asymm etric milharcsecon d rad io outflows (extend ing up 
. to -1000 AUl lMartf et all ( Il998l) and lParedes et all 12000 ) sug- 
' gested the presence of a mildly re l ativist ic (v ~ 0.2 c) jet. 
iParedes et atl (12002) and Ri bo et al.l (|2008|) have therefore pro- 
posed the source to be a microquasar. 

LS 5039 has been observed by severa l X-ray satellites. A 
summary of all observations can be found in lBosch-Ramon et al.l 
(2005). Flux variations on timescales of days as well as mini- 
flares on shorter timescales have been observed. RXTE observa- 
tions, covering the whole orbital period, show variability with a 
maximum of emission at orbital phase = 0.8, interpreted due to 
periastron passage. Moreover, the power-law shaped X-ray spec- 
trum appears harder in the high flux state (Bosch- Ramon et al.L 
I2005h . 

Because inhomogeneo us jets could emit y- rays via inverse 
Compton scattering, Pared eset al.l (l2000l |2002|) have proposed 



* Based on observations with INTEGRAL, an ESA project with in- 
struments and science data centre funded by ESA member states (espe- 
cially the PI countries: Denmark, France, Germany, Italy, Switzerland, 
Spain), Czech Republic and Poland, and with the participation of Russia 
and the USA. 



the source to be the counterpart of the unidentified hig h energy 
EGRET source 3EG J1824-1514 (Har tman et all 1 1999) and also 
a component of the MeV emitter GR O J 1823- 12, wh ich is a su- 
perposition of three EGRET sources (ICoUmaiil2003l) . 

The source has likewise been discovered to emit at TeV en- 
ergies with the H igh Energy Stereoscopic System (H.E.S.S.) 
dAharonian et al.L [2005). This demonstrated that X-ray bina- 
ries are capable of accelerating particles to TeV energies. The 
H.E.S.S. measurements have revealed that the flux and energy 
spectrum of the source are indeed mo dulated with the ~ 3.9 d 
orbital period of the binary system dAharonian et al 1 120061) . 
The VHE y-ray emission is largely confined to half of the or- 
bit and peaks around the inferior conjunction (INFC) epoch of 
the compact companion, where a hardening of the spectrum is 
observed. According to the authors these findings indicate that 
y-ray absorption with pair production occurs in the system. 
Interestingly, orbital variability at TeV y-ray energies has also 
been discovered in t he similar system LS I-h61°303 with MAGIC 
(lAlbertetaLL 120061) and VERITAS (.Acciari et al.Ll2008i) . 

The nature of the compact object in the system RX J1826.2- 
1450 /LS 5039 as well as the orgin of the TeV emis- 
sion is still under debate. A microquasar scenari o (see 
e.g. [Pared es et al., 2000; Bosch-Ramon & Paredei, 12004 
iTorres & Sierpows ka-Bartosik, 2008) as well as a pulsar sce- 
nario (see e.g iDubus, 2006; S ierpowska-Bartosik & Toiresl 
120071; iBosch-R amon et al., 2008) have been discussed. In this 
Letter, based mainly on archival data from INTEGRAL which 
span 4.6 years between March 2003 and October 2007, we 
confirm the detection of the source reported bv lBird et al.l d2007h 
in the energy range from 25 keV up to 200 keV. We report the 
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Table 1. Summary of selected and analyzed INTEGRAL data 
for LS 5039. Analyzed data were taken from 69 revolutions, i.e. 
from revolution 0050 (03-11-2003) to revolution 0614 (10-25- 
2007). 



z [R.] 

Fig. 1. A sketch of the system. The orbital values are listed in 
Casares et al.l (l2005 h and the definition of phase intervals is from 
Aharonian et all (12006) . 



discovery of the modulation of the hard X-ray flux with orbital 
phase. 



2. Observations and data analysis 

2.1. Observations 

The Internati onal Gamma-Ray Astrophysics Laboratory 
(INTEGRAL, IWinkler et akl l2003h carries on board two main 
gamma-ray inst ruments, the Spectrom eter on INTEGRAL (SPI, 
20 keV-8 MeV, IVedrenne et"alll2003l) . and the Ir nager on-Board 
INTEGRAL Satellite (IBIS, 15keV-10MeV, lUbertini et all 
I2003h as well as two monitoring instruments in the X-ray and 
optical rang e, the Joiri t Euro pean X-Ray Monitor (JEM-X, 
3-35 keV, Lund et al.', l2 003h. and the Optical Monitoring 
Camera (OMC, Mas-Hesse "etan l2003l) . 

1648 science windows (ScWs) in the region that in- 
cludes LS 5039 are found in the data archive avail- 
able at the INTEGRAL Scientific Data Center (ISDC, 
http://isdc.unige.ch/), including public data as well as data from 
the Key Programme. These ScWs were preselected choosing an 
off'-axis angle of the source smaller than 10 degrees and choos- 
ing only pointing observations that had a nominal observation 
time of more than 1000 seconds. 1627 ScWs were successfully 
processed through the analysis with OS A 7.0 up to the image 
step IMA. 

In order to reduce statistical effects a further selection of 
"good science windows" was necessary. ScWs are included in 
the analysis when they fulfill the following criteria: 

1. The effective exposure time after extraction is larger than 
1000 sec. 

2. The number of good time intervals (GTI) does not exceed 
10. 

3. The mean background intensity level of a single ScW in the 
energy range 25 - 60 keV falls within 3cr of the distribution 
of the mean background intensity levels of all other ScWs in 
the same energy range. 

The resulting 1393 ScWs cover the time range from 
03-11-2003 (ScW: 005000080010) to 10-25-2007 (ScW: 
061400950010). Firstly, the sample was split into two parts ac- 
cording to the orbital phase for imaging analysis, and secondly 



Number of all extracted ScWs 


1648 


Number of ScWs after selection 


1393 


Number of ScWs in phase SUPC 


721 


Number of ScWs in phase INFC 


672 


final total analyzed exposure time 


^ 2.99 Msec 


splitted in SUPC phase 


« 1.56 Msec 


and in INFC phase 


« 1.43 Msec 



into 10 smaller parts according to the orbital phase to analyse the 
flux dependence on orbital phase (see Sect l2.2l ). Table[T]summa- 
rizes the observational data of the source LS 5039 taken by the 
INTEGRAL satellite and the selection made for this analysis. 
Only in less than 10% of the selected ScWs was LS 5039 in 
the fully coded field of view of JEM-X. Based on the sensitivity 
of JEM-XQ and the amount of available data, a JEM-X analy- 
sis for a faint source like LS 5039 is not suitable. In this paper 
we therefore focus our analysis on the data of the Integral Soft 
Ganmia-Ray Imager (IBIS/ISGRI). 



2.2. Analysis 

For the analysis of ISGRI data the Off-line Science Analysis 
(OSA) software (version 7.0) provided by the ISDC was used. 
In addition, we cross checked the results and also generated the 
lightcurve in Sec. l2^.2lwith the analysis package developed by 
INAF-IFC PalermcPI (S egreto & Ferrigiiol 12006,) . 

A special handling of the OSA was used for this analysis. 
First, we used a user catalogue which includes 115 sources, 
i.e. all sources listed in the INTEGRAL General Reference 
Catalog which are in the field of view (FoV) and which turned 
out to be as bright as or brighter than RX J 1826.2- 1450 (= 
LS 5039). This is required because the source of interest is quite 
faint and the crowded field produces too many ghost images if 
the known sources are not subtracted. Second, the spectrum was 
produced with the help of the OSA tool mosaic_spec3 from 13 
mosaics, each covering one energy bin in the 13-250 keV range. 
The binning was chosen identical to the default energy binning 
OSA uses for the spectrum calculation step (SPE). This proce- 
dure is needed due to the fact that the source is faint and the 
FoV is too crowded, and provides more reliable results than the 
standard OSA pipeline does. 

We also checked our results using the standard SPE and the 
lightcurve creation (LCR) steps of OSA. In this case we used a 
special catalogue including all sources with a detection signif- 
icance greater than 6 cr in the energy band 20 - 60 keV. The 
catalogue includes 65 sources. According to the User Manual 
of the INTEGRAL IBIS/ISGRI data analysis this is beyond the 
standard method (see the IBIS Analysis Manual' for details). 



at: 
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' The User Manual to the OSA can be found 
http://isdc.unige.ch/?Support+documents See p. 15 
JEM-X Analysis User Manual for a JEM-X sensitivity plot. 

- http : //www ■ pa ■ ias f ■ cnr . it/~f errigno/ 1 
INTEGRALsof tware . html 

^ The INTEGRAL General Reference Catalogue is available at 
ISDC: http : / /isdc . unige . ch/index . cgi?Data+catalogs 

http : //isdc . unige . ch/Soft/download/osa/osa_doc/] 
osa_doc-5 . 8/osa_um_ibis- 5 . 8/nodel08 . html#mosaspec 
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However, our results show that a careful treatment of this stan- 
dard method provides a comparable result in spite of the restric- 
tions (see Fig.|4]and Sec. 12.2. 21 ). 

2.2.1. Imaging 

The source is clearly seen at a significance level of 7.7 cr in the 
lower energy band of 25 - 60keV in the significance map based 
on the sum of all ScWs fulfilling the selection criteria. In addi- 
tion, we divided the data set into two phase intervals: the superior 
conjunction passage at orbital phase 0.9 < (f> < 1.45 (SUPC), 
and the inferior conjunction, at orbital phase 0.45 < (f> < 0.9 
(I NFC) (see Figs.[T | and[4| for details). Phase intervals are defined 
in lAharonian et all ( |2006|) for the H.E.S.S. data. The two orbital 
resolved significance maps in the energy range 25 - 60keV are 
shown in Fig. ID Positions of sources (open boxes) are provided 
by the latest INTEGRAL General Reference Catalog^ 

At the INFC phase interval, the source was detected with 
a significance ^ 1.% cr, for an exposure time of ^ 1.43 Msec 
(Fig. |2] lower panel). At the SUPC phase interval the source 
was observed for a total exposure time of 1 .56 Msec, and was 
only marginally detected: at the source position there is an ex- 
cess at the level of 2.8 cr (Fig. |2l upper panel). For the total 
sample and the total energy band 25 - 200 keV the source is 
detected at the position fitted by OSA (RA - 276.517 deg, 
DEC - -14.812 deg) with 7.4cr. Within the e rrors, this is consis- 
tent with the source position determinated bv lBird et al.l (l2007h . 

2.2.2. Timing analysis 

A lightcurve pro duced w ith the INAF-IFC data analysis soft- 
ware ( Seg reto & Fe rrignd 12006) was used for the timing analy- 
sis. According to the authors this software is suitable to produce 
reliable results even in the case of the presence of a larger num- 
ber of sources in the FoV. 

Searching for an unambiguous period failed due to low 
statistics and large gaps in the sample of observations. But a pe- 
riodogram produced by epoch folding provides a peak at a pe- 
riod of 3.903 days, which is within uncertai nties i dentical to the 
orbital period determined by Casares et al. (|2005|) in the optical 
band. Fig.|3]shows the periodogram. 

In a next step the lightcurve was folded using the known 
orbital period of Porb = 3.90603 days with HJD (Tq) = 
2451 943'?09 ± 0?10 (= phase ^ 0.0 and periastron passage) 
determined bv lCasares et al.l (l2005l) . to produce a 10 bin orbital 
phase profile. In addition, 10 mosaics corresponding to each of 
the 10 orbital phase intervals were produced and the count rates 
at the catalogue^ position were determined. The values of the 
fluxes were derived from a 2D-Gaussian fit with fixed width 
(HWHM - 6.0 arcmin) and fixed source position. Both methods 
yield (within errors) comparable results: The hard X-ray emis- 
sion is clearly modulated with the orbit, and appears to be in 
phase with the TeV emission (Fig.j^. The H.E.S.S. data points 
in the figure are from lAIiaronian et al.l ( 120061) . 

2.2.3. Spectral analysis 

A statisticaUy significant IBIS/ISGRI spectrum of LS 5039 was 
obtained for the orbital phase interval corresponding to INFC 
(not shown here) and was modeled with a power law with F - 
2.0!° 2 (90% confidence level). The reduced for the fit is 1 .4 (9 
dof). Mean countrates for the two phase intervals allow us to es- 
timate a flux for the INFC phase interval of (3 .54 + 2 .30) x 1 0" " 
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Fig. 2. INTEGRAL significance maps in the energy range 25 - 
60 keV. The open boxes are known source positions. Comparable 
exposure times and a higher significance for LS 5039 in the 
INFC phase interval indicates variability of the hard X-ray emis- 
sion with orbital phase (see also Fig.|4|. The INFC phase inter- 
val includes the maximum flux state of the source. Details can 
be found in the text. 



erg cm"^ s ' (90 % conf. level). For the SUPC phase interval a 
flux upper limit of 1.45 X 10"'^ erg cm~^ s"' (90 % conf. level) 
for the energy range 25 - 200 keV can be determined. 



3. Summary and discussion 

In this paper we report on a 3 Msec INTEGRAL monitoring 
of the HMXRB RX J 1826- 1450/ LS 5039. The faintness of the 
source and the fact that it lies in a crowded field required a very 
careful, non-standard handling of the data. 

Our results, obtained by applying diff'erent analysis methods, 
point out that the source significantly emits at hard X-rays (25 - 
200 keV) and that this emission varies with the orbit in phase 
with the very high energy y-rays detected with H.E.S.S. (Fig. 
|4]l. The spectrum at the inferior conjunction is well described 
by a power law, while at the superior conjunction the hard X- 
ray emission, if any, is below the sensitivity of INTEGRAL. Our 
result indicates that accretion might not be the mechanism for 
the production of the hard emission, since, in this case, we would 
expect a rather sharp flux maximum near periastron. 

Moreover, the observed close correlation of hard X-ray and 
VHE y-ray emission in phase suggests that the emission would 
originate from the same source region, possibly from particles 
produced in the same acceleration process. 
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Fig. 3. An epoch folding on the INTEGRAL lightcurves shows 
a significant peak at a period that is nearly identical to the or- 
bital period known and determi ned from optical observation of 
the system dCasares et all 12005). All other peaks are likely beat 
frequencies due to the large gaps in the sample of observations. 



This casts doubt on the simple explanation that the TeV 
variability is mainly caused by y-y absorption, as suggested in 
IXharonian et al. (2006). 

We note that, from the correlation alone, more compli- 
cated scenar ios involving y-y absorption such as described by 
iDubus et~aLl (^008) cannot be ruled out. However the cascades 
initiat ed in the y -y absorption process produce secondary elec- 
trons dBednareM, [l997) also emitting X-ray synchrotron emis- 
sion. Therefore, soft and hard X-ray flux levels can be used to 
constrain the amount of absorption involved. The detection of 
TeV emission at SUPC together with the upper limit in the hard 
X-ray domain derived here may pose problems for models in- 
voking only one emi ssion region close to t he com pact object, as 
recently discussed in lBosch-Ramon et al.l (|2008[) . 

Radio observations interpreted as jets have led to a micro- 
quasar interpretation for this source. In this scenario, particles 
are accelerated to k TeV energies jPerucho & Bosch-Ramonl 
2008i) due to the interaction bet ween the jet an d the stellar wind 
of the companion. Alternativelv. lDubusI (l2006l) has suggested the 
source to be a pulsar binary system, where VHE emission arises 
from the pulsar wind material shocked by the interaction with the 
stellar companion wind and afterwards flowing away in a comet- 
shape tail. Unfortunately, the observations presented here do not 
enable us to discriminate betwee n the two scenarios , although 
the predicted X-ray variability bv lParedes et alj ( 120061) seems to 
be not in accordance wit h our observations. Finally, we remark 
that a hadronic scenario dRomero et all 120031; iRomero & Vilal 
|2008|) is not ruled out by our findings. A longer and more de- 
tailed interpretation of our results is beyond the scope of this 
Letter, and will be the subject of further investigations. 
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Fig. 4. Top: Orbital profile for LS 5039 in the energy range 
25 - 200 keV observed with INTEGRAL. The lightcurve is 
folded with the orbital period of 3.9 days. Bottom: Orbital pro- 
file for L S 5039 for the energy ra nge above 1 TeV observed with 
H.E.S.S. dAharonian et al I I2OO6I) . The hard X-ray emission is in 
phase with the TeV y-rays. 
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